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In this work, the ab initio calculations have been carried out to study the transition proba- 
bilities corresponding to the forbidden transitions of astrophysically important electromag- 
netic transitions in the singly ionized zinc, Zn II. Many important electron correlations are 
considered to all-orders using the relativistic coupled-cluster theory. Calculated ionization 
potentials are compared with the experimental values, wherever available. To our knowledge, 
oscillator strengths of the magnetic dipole and electric quardupole transitions are estimated 
for the first time. The transitions span in the range of ultraviolet, visible, and near infrared 
regions and are crucial for the astrophysical observations. 

I. INTRODUCTION 

With the advancement of satellite based observations, we can get the information about the 
abundance of elements in the interstellar space which provides a probe of the nature and composi- 
tion of the interstellar medium. The data obtained with the Goddard High Resolution Spectrograph 
(GHRS), the Hubble Space Telescope (HST) and other satellites based on the spectroscopic facili- 
ties exhibit very high signal-to-noise ratios and allow routine measurements of very weak lines (e.g. 
forbidden lines) with attainable accuracy. 

Generally allowed electric dipole transitions (El) are important because of their strengths 
and easily accesible for the incoming observations. However, forbidden transitions, like electric 
quadrupole (E2) and magnetic dipole transitions (Ml), have also great importance in the astro- 
physics even though they are weak, because they carry valuable informations like thermal doppler 
effects of heavenly bodies. The forbidden lines provide important clues in other areas of astro- 
physics, beacuse of the long lifetime of the upper state against radiative decay. These lines are 
particularly sensitive to the collisional de-excitation and serve as indicators of electron density and 
temparature, and Tg, in the emmision region. Determination of and Tg from the forbidden 
line intensities was discussed origanally for the general case by Seaton [l| and Seaton and Oster- 



brock [2]. A number of such transitions have been observed in the ultraviolet spectrum of the solar 



y observed in quasars with an intensity 
^. Moreover, gaseous nebulae exhibit in 



corona. Forbidden atomic emmision lines are common ; 
often comparable to accompanying 'allowed transition' 
their spectra forbidden transition lines of low excitation energy. Many astrophysical phenomena 
like coronal heating, evolution of chemical composition in stellar envelopes, determination of the 
chemistry in the planetary nebulae precursor's envelope are believed to be explained largely by 
these forbidden lines. In laboratory tokamak plasmas and in various astronomical objects, suit- 
ably chosen these forbidden lines serve as a basis for reliable electron density and temperature 
diagnostics 4 1 . 

Zn has been idetified in Be stars [5|] and experimental data for absorption lines from the ground 
level are also available. Sneden et al. had studied the abundance of Zinc (Zn) as a function 
of metallicity for a number of metal poor stars. Zn is closely tracking ion in abundance for most 
stars. There are peculiar stars for which Zn is either scarce (if not non-existing) or over abundant 
[3]. An UV-atlas with line identifications [sl is now available for the HgMn star x-Lupi jol, which 
has been studied extensively. The signature of singly ionized zinc (Zn II) are almost undetectable 
in these stars, apart from the interstellar medium (ISM). Zn II relatively undepleted in the diffuse 
interstellar medium, and hence may be used to infer the amount of H in sightline. The lines of Zn II 
have been observed in the interstellar absorption spectrum towards C Ophiuch i [lo| . The analysis 



111 ] in determining 



data of Zn II spectrum from GHRS on the HST are done by Savage et al. 
the depletion of Zn II. Allowed transitions and corresponding oscillator strengths of this ion have 
been studied by atomic physicsts, both experimentalists and theorists Ij], but there has been no 



work towards the 'forbidden transition' of Zn II till date in our best of knowledge, which are very 
important in astronomy. 

Here, we have employed the relativistic coupled cluster method with single and double excita- 
tions (RCCSD) for determinating the 'forbidden transition' amplitudes. This is one of the most 
powerful highly correlated approaches due to its all order behavior for the correlation operator 
13l |. The Fock-space multi-reference coupled cluster (FSMRCC) theory for one electron attach- 



16l |. We provide a brief review of 



ment process used here has been described elsewhere 
this method. The theory for a single valence system is based on the concept of common vacuum for 
both the closed shell N and open shell Nil electron systems, which allows us to formulate a direct 
method for energy differences. Also, the holes and particles are defined with respect to the common 
vacuum for both the electron systems. Model space of a (n,m) Fock-space contains determinants 
with n holes and m particles distributed within a set of what are termed as active orbitals. For 



3 



example, in this present article, we are dealing with (0,1) Fock-space which is a complete model 
space (CMS) by construction and is given by 



i 

where Cj^'s are the coefficients of ^^"'^^ and (^■'^'^^'s are the model space configurations. The 
dynamical electron correlation effects are introduced through the valence-universal wave-operator 

n = {exp(5)} (1.2) 

where 

m n 

S = J2J1 S'^"''^ = S^''""^ + S^'''^ + -S^''"^ + • ■ • (1-3) 

A:=0 ;=0 

At this juncture, it is convenient to single out the core-cluster amplitudes S^^'^^ and call them T. 
The rest of the cluster amplitudes will henceforth be called S. Since is in normal order, we can 
rewrite Eq. ()1.2p as 

Q. = exp{T){exp{S)} (1.4) 

In this work, single (Ti,5i) and double excitations (T2,5'2) are considered for T and S clusters 
operator. Wavefunction of the system with single valence orbital v, therefore, yields the form 

= Q.\^df) = e^i+^Hl + Siv + S2v}\^df). (1.5) 

Triple excitations are included in the open shell RCCSD amplitude which correspond to the corre- 
lation to the valence orbitals, by an approximation that is similar in spirit to CCSD(T) 0. The 
approximate valence triple excitation amplitudes are given by 



abk — . . . . ^ _ . ) 



abk — ' (1'^) 

~t~ ~t~ £k 

where S^^'^'^^H. are the amplitudes corresponding to the simultaneous excitations of orbitals a, b, k 
to p,q,r, respectively; VT2 and VS^^'^\ are the connected composites involving V and T, and 
V and S^^'^\ respectively, where V is the two electron Coulomb integral and e's are the orbital 
energies. 
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The the transition element due to any operator D can be expressed, in the RCCSD method, as 

{^f\D\^i) 



D 



'^f\{l + Sf^e''^DeT{l + Si}\1>i) 



The one-electron reduced matrix elements of Ml and E2 operators are given by 



(1.7) 



1^. 



6 Ki + Kf 

ok 2 



drji{kr){Pf{r)Q,{r) + Qf{r)Pi{r)) 



(1.8) 



and 



15 



X 



drj2{kr){Pf{r)Pi{r) + Qf{r)Qi{r)) + h{kr){ 



+ {Pf{r)Qi{r)-Qf{r)n{r)) 



K f Ki 



){Pf{r)Q,{r) + Qj{r)n{r)) 
(1.9) 



respectively. Here, ji and Ki = i: (ji + are the total angular momentum and relativistic angular 
momentum quantum numbers, respectively, of the 'f^ electron orbital. The quantity Cm is the 
Racah tensor and ji{kr) is the spherical Bessel function of order /. P^j and Q/^i are the large and 
small radial components of the Dirac wavefunctions. 

The emission transition probabilities (in sec ^) for the E2 and Ml channels from states / to z are 
given by 



^ 1.11995 xlQig ^^, 



A 



Ml 



X^2jf + 1 
2.69735 X 10^3 



S 



Ml 



[1.10) 



(1.11) 



''^ A3(2j/ + 1) 

where S = |(^'/|-D|^'j)p is the transition strength for the operator D (in a.u.) and A (in(A)) is the 
corresponding transition wavelength. 



II. RESULTS AND DISCUSSIONS 



We calculate the DF wavefunctions |$df) using the Gaussian- type orbitals (GTO) as given in 



25l | using the basis functions of the form 



= Cj/Ve-^^^' (2.1) 
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with A; = 0, 1,2,3, .... for s, p, d, f type orbital symmetries respectively. The radial functions 

'F^' and ^F^' represent the basis functions correspond to large and small components of the Dirac 

L /s 

orbitals. are the normalization constants which depend on the exponents. The universal 

even tempering condition has been applied to the exponents ; i.e., for each symmetry exponents 
are assigned as 

ai = ao/3*"^ « = 1,2, (2.2) 

where N is the number of basis functions for the specific symmetry. In this calculation, we have 
used oq = 0.00831 and /3 = 2.99. The number of basis functions used in the present calculation is 
32, 32, 30, 25, 20 for / = 0, 1, 2, 3, 4 symmetries, respectively. 

Number of DF orbitals for different symmetries used in the RCCSD calculations are based on 
convergent criteria of core correlation energy for which it satisfies numerical completeness. There 
are only 11, 10, 10, 7 and 5 active orbitals including all core electrons are considered in the CCSD(T) 
calculations for / = 0, 1, 2, 3, 4 symmetries, respectively. We first calculate T amplitudes using the 
RCCSD equations of closed shell system (Zn II) and then solve the S amplitudes from the open 
shell equation for this single valence states of Zn II. 

In Table I, we present the ionisation energy and fine structure splittings and their comparision 



with the NIST 26|] results. Our calculated ionisation energies are in excellent agreement with 
the NIST values. The present study uses Dirac Hartee-Fock (DHF) orbitals to form the slater 
determinants to represent different atomic states and includes the effects of Coulomb correlation 
using a fully ab initio all-order many-body RCCSD method. 

The oscillator strengths for transitions between different energy states depend linearly on energy, 
and quadratically on the matrix elements as shown in Eq. (7). This shows the necessity to calculate 
transition matrix elements appropriately. 

In table III, we have given the ab initio transition probabilities {Afi value) for the 'ground 
to excited' and 'excited to excited' states transitions, which are also astrophysically important. 
Strong transition probabilities have been observed for many E2 transitions. 



"Laser-produced plasma [12 
^LIF [27] 

"LIF -|- Hollow cathode lamp 



'^Phase shift 
'^Beam foil 



29] 



28| 



•^Relativistic many-body perturbation 3l|] 
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TABLE I: Ionization Potentials (IPs) and fine structure (FS) splittings (in cm^^) of Zn II and their compar- 
ison with NIST values. The percentage of differences of our calculated IP results compared to NIST results 
are shown in parenthesis on the side of CCSD(T) IP results. 

IP FS 
States CCSD(T) NIST CCSD(T) NIST 



4si/2 


000 


00 


000 


00 


000.00 000.00 


4pi/2 


48557 


23 


48481 


00 




4p3/2 


49422 


75 


49355 


04 


865.52 874.04 


5Sl/2 


88379 


53 


88437 


15 




4^3/2 


96781 


13 


96909 


74 




4^5/2 


96830 


94 


96960 


40 


49.81 50.66 


5pi/2 


101329 


53 


101365 


93 




5P3/2 


101623 


61 


101611 


43 


294.08 245.50 


6si/2 


114617 


73 


114498 


02 




4/5/2 


117026 


80 


117263 


40 




4/7/2 


117026 


39 


117264 


00 


0.41 0.60 


5^3/2 


118112 


35 


117969 


32 




5^5/2 


118140 


33 


117993 


61 


27.98 24.29 


6Pl/2 


119848 


71 


119888 


51 




6P3/2 


119944 


98 


119959 


34 


96.27 70.83 



3MCDF + Core polarizationjsal 
'^RHF + Core polarizatior i [3^ 
* Coulomb approximation [3j] 

Relativistic Supersymmetry Quantum Defect Theory 
'^One electron approx. with core polarisation [3r' 
'WKB (sinusoidaloscillation) 
"^CI + Core polarization 



37] 



"One electron approx model potential 
°HF [401 



PNIST [261] 
■^Beam foil 
''Beam foil 
*Beam foil 



41] 
42] 
43] 
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TABLE II: Oscillator strength corresponding to allowed transitions among different low-lying states of Zn 
II. 



State 




Experiment 


Other theories 


This work 


4si/2- 


^ 4pi/2 


0.249(50),"^ 0.249(29),^ 0.256(26),^= 
0.209(18), 0.306(58), 0.246,^ 


0.317," 0.291,'' 0.290,^= 
0.260,-' 0.2526,3 0.238," 

0.2965,* 0.2521,J 0.2,^= 
0.2856, 0.262,™ 0.268, 


0.2621 


4si/2- 


* 4p3/2 


0.467(93)," 0.48(55),'' 0.488(49),^= 
0.513(60),^ 0.606(53), 0.62(6) 
0.590(84),^= 0.406(41)," 0.501,^ 


0.642," 0.586,'' 0.615,*^ 
0.515,-' 0.5159,3 0.487,"^ 
0.6046,' 0.5187,^ 0.476,° 
0.41, 0.5o21, 0.53/, 0.552, 


0.535 


4pi/2- 


44/2 


0.543(^0), 0.000,'^ 


A AC1 a A 70/iOA 

0.951, 0.78489,'' 

A lOCi h A /I ^'y i 

0. /89, 0.84/ /, 


A 0010 

0.8218 


4p3/2- 


^ 44/2 


U.Uo I (11), 


A AOQC; d A A7ACO Q 
A A'yAA h A AOCAO i 

U.0/99, 0.08593, 


U.UoDo 


4p3/2- 


^ 44/2 


0.69(14)," 0.684(57),'' 
0.704(75),'' 0.258(47),'' 0.555,?' 


0.843," 0.757,'' 0.71483,-'' 
0.718,'' 0.7721,' 0.67,° 


0.7823 


44/2- 


^ 4/5/2 


0.99(20)," 0.868,P 


1.08," 1.057,* 


1.086 


44/2- 


^ 4/7/2 


1.32(26)," 


1.02," 1.009,* 


1.131 


5pi/2- 


^ 54/2 


0.82(0.16)," 


0.826," 0.9428,* 


0.9541 


5p3/2- 


■> 54/2 


0.0240(48)," 


0.0787," 0.09629,* 


0.0949 


5^3/2- 


■> 54/2 


0.52(10)," 


0.708," 0.8633,* 


0.5706 



III. CONCLUSION 

The continuing developments in astrophysical and astronomical observations demand accurate 
theoretical transitions data to determine stellar chemical composition. Highly correlated relativistic 
coupled cluster theory has been employed to study the oscillator strengths of the astrophysically 
important forbidden transitions, Magnetic dipole and electric quadrupole transition amplitudes 
among the bound states of Zn II are important for astronomical observations and plasma researches. 
To our knowledge these are the first calculations of oscillator strengths for the forbidden transitions 
of Zn II. All the transitions are in the ultraviolet, visible or near infrared regions. This work will 
motivate astronomers to observe these lines of Zn II to predict the abundances of these species in 
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astronomical bodies and experimentalists to verify our results. 
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TABLE III: Transition wavelengths and transition probabilities of Zn II. 



Terms Acc5d(A) A/j-value E2 ^/j-value Ml 



44/2 


^ 44/2 


2007588.68 


2.3078E-12 


1.326E-06 




54/2 


4687.96 


4.7206E+01 


5.6123E-01 




^ 54/2 


4681.82 


1.4785E+01 


3.4344E-04 




^ 4Si/2 


1031.88 


2.1426E+04 






^ 5Si/2 


11902.49 


2.4488 






6Si/2 


5606.44 


3.2714E+01 




44/2 


^ 54/2 


4698.93 


2.2248E+01 


4.0621E-04 




54/2 


4692.76 


5.3873E+01 


2.9588 




^ 4si/2 

/ 


1032.72 


3.1598E+04 






^ 5Si/2 


11832.34 


3.7945 






^ 6Si/2 


5612.14 


4.9539E+01 




54/2 


^ 54/2 


3573623.79 


2.4516E-12 


2.3577E-07 




4si/2 


846.65 


4.4728E+03 






^ 5Si/2 


3363.28 


6.8785E+02 






^ 6Si/2 


28615.41 


5.1564E-01 




54/2 


^ 4si/2 


846.45 


6.5596E+03 






^ 5Si/2 


3360.12 


1.0266E+03 






^ 6Si/2 


28388.95 


8.0749E-01 




4pi/2 


4p3/2 


115537.36 


1.3904E-06 


5.8297E-03 




^ 5pi/2 


1894.93 




3.6309 




^ 5p3/2 


1884.43 


5.6262E+02 


9.8766E-02 




6pi/2 


1402.69 




1.6186 




^ 6p3/2 


1400.80 


1.5624E+02 


5.4193E-02 




4/5/2 


1460.51 


3.3301E+03 




4p3/2 


^ 5pi/2 


1926.53 


1.1498E+03 


2.8073E-01 




5p3/2 


1915.67 


4.8138E+02 


6.4942E+01 




^ 6pi/2 


1419.93 


2.9788E+02 


9.9683E-02 




^ 6p3/2 


1417.99 


1.1143E+02 


3.1170E+01 




^ 4/5/2 


1479.20 


9.6132E+02 






^ 4/7/2 


1479.21 


4.3255E+03 




5pi/2 


5p3/2 


340039.01 


1.8059E-07 


2.2860E-04 




^ 6pi/2 


5399.81 




1.1853E-01 




^ 6p3/2 


5371.88 


5.6369E+01 


2.3181E-03 




^ 4/5/2 


6370.53 


9.7419E+01 
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TABLE IV: Continuation from table III 



Terms Accsd(A) A/j-value E2 ^/j-value Ml 

5p3/2 ^6pi/2 5486.96 1.1319E+02 1.5547E-02 

2.2892 



6pi/2 
6p3/2 
4Sl/2 
5Sl/2 



6pi/2 
6p3/2 

4/5/2 
4/7/2 

6p3/2 
4/5/2 
4/5/2 
4/7/2 
5Sl/2 



6s 



1/2 



6Sl/2 



5486.96 1.1319E+02 

5448.11 4.6981E+01 

6492.16 2.5774E+01 

6492.33 1.1596E+02 
1038663.19 7.3181E-09 
35437.11 2.8051E-02 
34267.93 9.3366E-03 
34263.15 4.2045E-02 
1131.48 
872.46 
3811.23 



8.0184E-06 



1.2592E+02 
4.9760E+01 
2.6961 



